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The creation of a quantum Universe is described by a density matrix which yields an ensemble 
of universes with the cosmological constant limited to a bounded range Amin < A < Amax- The 
domain A < Amin is ruled out by a cosmological bootstrap requirement (the self-consistent back 
reaction of hot matter). The upper cutoff results from the quantum effects of vacuum energy and 
the conformal anomaly mediated by a special ghost-avoidance renormalization. The cutoff Amax 
establishes a new quantum scale - the accumulation point of an infinite sequence of garland-type 
instantons. The dependence of the cosmological constant range on particle phenomenology suggests 
a possible dynamical selection mechanism for the landscape of string vacua. 

PACS numbers; 



Quantum cosmology [l|, 0| and Euclidean quantum 
gravity Q might effectively restrict the landscape of 
string vacua. This landscape is too big Q to predict ei- 
ther the observed particle phenomenology or large-scale 
structure formation within string theory itself. Other 
methods have to be invoked, at least some of them based 
on the cosmological wavefunction 0, 0, • This is domi- 
nated by the exponentiated Euclidean action, exp(— S'e), 
calculated on the gravitational instanton which is a sad- 
dle point of an underlying path integral over Euclidean 
4-geometries. This instanton gives rise to Lorentzian sig- 
nature spacetime by analytic continuation across mini- 
mal hypersurfaces. The continuation can be interpreted 
either as quantum tunneling or as the creation of the Uni- 
verse from "nothing" . Thus, the most probable vacua of 
the landscape become weighted by the minima of S'e- 
This might serve as a method of selecting a vacuum from 
the enormously big string landscape. 

An immediate difficulty with this program arises from 
the infrared catastrophe of small cosmological constant 
A. The Hartle-Hawking wave function 3], which de- 
scribes nucleation of the de Sitter Universe from the Eu- 
clidean 4-dimensional hemisphere, has the form 



^-HH ~ exp(-S'E) = exp(37r/2GA). 



(1) 



This diverges for A — > because of unboundedness of 
the Euclidean gravitational action. Despite some early 
attempts to interpret it as the origin of a zero value of 
A Q, this result remains both controversial and anti- 
intuitive because it disfavors inflation and prefers cre- 
ation of infinitely large universes. Apart from the tun- 
neling proposals of 9] which employ an opposite sign in 
the exponential of ([T]) and thus open the possibility for 



opposite conclusions [lO| , no convincing resolution of this 
problem has thus far been suggested. 

In this Letter we show that Euclidean path integration 
framework naturally avoids this infrared catastrophe. We 
attain this result by: i) extending the notion of Hartle- 
Hawking pure state to a density matrix which describes 
a mixed quantum state of the Universe and ii) incorpo- 
rating the nonperturbative back reaction of hot quantum 
matter on the instanton background [ll|. These exten- 
sions seem natural because whether the initial state of 
the Universe is pure or mixed is a dynamical question 
rather than a postulate. We address this question below 
by embedding both types of states into a unified frame- 
work of a density matrix. 
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FIG. 1: Density matrix instanton. Dashed lines depict the 
Lorentzian Universe nucleating at minimal surfaces E and E'. 



A density matrix p[ tp, ip'] is represented in Euclidean 
quantum gravity [l^ by an instanton having two disjoint 
boundaries S and S' associated with its two entries and 
tp' (collecting both gravity and matter observables). The 
instanton interpolates between these, thus establishing 
mixing correlations, see Fig[T] In contrast, for the den- 
sity matrix of the pure Hartle-Hawking state the bridge 
between E and T,' is broken, so that the instanton is a 
union of two disjoint hemispheres which smoothly close 
up at their poles (Figl2|) — a picture illustrating the fac- 
torization of /5 = |4'hh)(^'hh|- 

The main effect that we advocate here is that thermal 
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FIG. 2: Density matrix of the pure Hartle-Hawking state rep- 
resented by the union of two vacuum instantons. 



fluctuations and quantum conformal anomaly destroy the 
Hartle-Hawking instanton and replace it with one filled 
by radiation. This is already manifest in classical theory, 
specifically in the Euclidean Friedmann equation for a 
scale factor a{T), a^ja^ = \/ — H"^ — C / (we use spa- 
tially closed FRW metric ds^ = N^{t) dr^ + a2(r) d^l^'^) 
in the gauge N = \ and express A = iH'^ in terms of the 
Hubble constant). The radiation density C/a^ prevents 
the half-instantons from closing and allows a to vary be- 
tween two turning points [H, a± = {^/ V^H){1 ± 
Vl — 4Ci?^)^/^. This forces a tubular structure on the 
instanton which spans (at least) one period of oscillation 
between a±, provided the constant C characterizing the 
amount of radiation satisfies the bound AH^C < 1. 

The existence of radiation naturally follows from the 
partition function of this state associated with the 
toroidal instanton obtained by identifying S' and E. The 
radiation back reaction supports the instanton geometry 
in which it exists. Remarkably, when the vacuum energy 
and conformal anomaly are taken into account this boot- 
strap yields a set of instantons - a landscape - only in 
the bounded range of A, 



A„ 



< A < Ai, 



(2) 



All values A < Amin are completely eliminated either 
because of the absence of instanton solutions or because 
of their infinitely large positive action. A similar situation 
holds for A > A^ax ~ no instantons exist there, and the 
Lorentzian configurations in this overbarrier domain (if 
any) are exponentially suppressed relative to those of ([2]). 

To quantify the above picture consider the density ma- 
trix given by the Euclidean path integral 



,r 



exp(-5E[.g,<^]), (3) 



where S-E[g,(t>] is the classical action, and the integration 
runs over gravitational g and matter fields interpolat- 
ing between ip and ip' at E and E'. The statistical sum 
exp{~r) is given by a similar path integral over periodic 
fields on the torus with identified boundaries S and E'. 

The back reaction follows from decomposing [g, 0] into 
a minisuperspace go{T) = (a(T), Af(T)) , and the "mat- 
ter" sector which includes also inhomogeneous metric 
perturbations on minisuperspace background (l>{x) ~ 



{(t){x), 'ilj{x),Af, (x), h^^{x), 



With a relevant decom- 



integral for F expresses in terms of the effective action 
r[gQ{T)] of quantized matter on the background goir), 
exp{-r[go]) = J D<P{x) exp(-5'E[5o, ^'(a;)]), as 



= j Dgoir) exp(-r[go(r)]). 



(4) 



Our approximation will be to consider r' [go (■'')] in the 
one-loop order, r[ga] = S-Eigo] + A-ioop[5o], and handle 
(jlj at the tree level, which is equivalent to solving the 
effective equations for -r[go]- 

Remarkably, r[gQ] is exactly calculable for confor- 
mally-invariant fields by a conformal transformation [l5l | 
relating generic FRW metric ds^ = a^ds^ to that of a 
static universe of a unit size, ds^ — dif' -I- d^Vl'^^^ (these 
metrics are denoted below as g and 5, while 77 is the con- 
formal time). The total action reads 

r[a{T),N{T)] =2 j ^ dT^-'^-Na + NH^a^^ 
+B f ^ dTN/a + f(2 f * drN/ay (5) 



We work in units of mp = ^37r/4G, and the integration 
runs between two turning points at t±. The first line is 
the classical part, the second line is the conformal con- 
tribution and the last line is the one-loop action on the 
static instanton of the metric g. 

The conformal contribution /^i_ioop[5] — A-ioop[5] is 
determined by the coefficients of OR, the Gauss-Bonnet 
invariant E ~ R^^^^^^—AE?^^+E? and Weyl tensor term in 
the conformal anomaly g^n Sri_ioop/Sg^.=g^/^{aDR + 
f3E + 7C^^^^)/4(47r)^. Specifically this contribution can 
be obtained by the technique of [l^; it contains higher- 
derivative terms ~ which produce ghost instabilities 
in solutions of effective equations. However, such terms 
are proportional to the coefficient a which can be put to 
zero by adding the following finite local counterterm 



^^[.9] = A-ioop[.g] 



2(47r)2 12 



d^xg'/'R^g). (6) 



position of the measure D[g,(j)] = Dgolr) x D<l>{x), the 



This ghost- avoidance renormalization is justified by the 
requirement of consistency of the theory at the quantum 
level. The contribution / i?[(7] — A? [5] to the renormalized 
action then gives the second line of ([5]) with B = 3/3/4. 

The static instanton with a period 770 playing the 
role of inverse temperature contributes -Ti-iooptff] = 
EoTjo + _F(ryo), where the vacuum energy Eq and free 
energy F{rio) are the typical boson and fermion sums 
over field oscillators with energies w on a unit 3-sphere 
^^0 = ±Eu^/2^ F{7^o) = ±j:jn{lTe''^^°). The 
renormalization ([6]) which should be applied also to 
A-ioop[5] modifies Eq, so that = Cq r/o + F{r]o), 

Co = Eq + 3a/16. This gives the third line of Eq.® 
with Co ~ B/2. This universal relation between Cq and 
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B = 3/3/4 follows from the known anomaly coefficients 
[l3| and the UV-renormalized Casimir energy in a static 
universe [l^ for scalar, Weyl spinor and vector fields re- 
spectively having: 
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17. (7) 



It is important that for conformally invariant fields the 
nonlocal action ([5|) is exact, and contains no other terms 
of higher order in the curvature. 

The effective equation 5r/5N{T) = has the form of 
the classical equation modified by the quantum i?-terni 



a 



c = B/2 + F'{m), m 



f/2 ^ 



+ 



dr / air) 



(8) 
(9) 



Remarkably, the contribution of the nonlocal Firjo) in ([5]) 
reduces to the radiation constant C as a nonlocal func- 
tional of a(r), determined by the bootstrap equation ([9]). 
Here F'{rio) = dF [rji^) / drjQ > is the energy of a hot gas 
of particles, which adds to their vacuum energy B /2. 




FIG. 3: Instanton domain in the (i/'^, C)-plane. Garland fam- 
ilies are shown for fe = 1, 2, 3, 4. Their sequence accumulates 
at the critical point (1/2B, B/2). 

Periodic instanton solutions of Eqs.fUl)-® exist only 
inside the curvilinear wedge of (i/2^ C')-plane between 
bold segments of the upper hyperbolic boundary and the 
lower straight line boundary of Fig|31 



4CH^<1, C > B - B^H^, BH^< 1/2. 



(10) 



Below this domain the solutions are either complex and 
aperiodic or suppressed by infinite positive Euclidean ac- 
tion. Above this domain only Lorentzian (overbarrier) 
configurations exist, but they are again exponentially 
damped relative to instantons in PH)) . 

These properties are based on the fact that the turning 
points of ([5]) exactly coincide with classical a-t, but a_ 



exists only when a^_> B, which gives rise to (|10p . Oth- 
erwise, a(r) at the contraction phase becomes complex 
or runs to zero which violates instanton periodicity. In 
the latter case a smooth Hartle-Hawking instanton with 
a_ = forms and yields 770 — *■ cxd in view of ([9]), so that 
F{vo) ~ P'ivo) ^ 0. Therefore, its on-shell action 



ro^Fi7jo)~VoF'ivo)- 



dad ( „ 



Bd' 



(11) 



due to B > diverges to -l-oo at a_ = and completely 
rules out pure-state instantons [Tl| . 

Moreover, inside the range PH)) our bootstrap elimi- 
nates the infrared catastrophe of A — > 0. Indeed rjo ~^ 00 
as 0, so that due to ^ C ^ B/2, but this is im- 

possible because in view of (fTO|) C > S at = 0. Thus, 
instanton family never hits the C-axes of — and 
can only interpolate between the points on the bound- 
aries of the domain pI7|) . For a conformal scalar field 
the numerical analysis gives such a family [11] starting 
at ff2 2.00, C !v 0.004, Fq « -0.16, and terminating 
at i?2 ^ 13 0, C « 0.02, Fo « -0.09. The upper point 
describes the static universe filled by a hot radiation with 
the temperature T ~ H/tt^/1 — 2BH'^, whereas the lower 
point establishes the lower bound of the A-range. 

The upper bound of the landscape follows from the 
existence of garlands that can be obtained by glueing 
tog ether into a torus k copies of a simple instanton [iSj, 
Il9j : see Fig|4l Their formalism is the same as above 
except that the conformal time in ^ and the integral 
term of (fTT|) should be multiplied by k. As in the case of 




FIG. 4: The garland segment consisting of three folds of a 
simple instanton glued at surfaces of a maximal scale factor. 

k = I, garland families interpolate between the lower and 
upper boundaries of (fTO|) . They exist for all fc, 1 < fc < 
00, and their infinite sequence accumulates at the critical 
point C = B/2, ij2 — 1/2B, where these boundaries 
merge. Within the l/A;2-accuracy the upper and lower 
points of each family coincide and read 
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-B- 



W fc2 

4fc27r2 ■ 



(12) 



With a growing k, garlands become more and more static 
and cold with T(fc) ~ l/(^/Blnfc2) _^ 0. Contrary to ^ 
the garland action is not additive in k, so that as fc ^ 00 
garlands do not dominate the ensemble. Their sequence 
converges to the instanton with -ff^^x = 1/2S, which 
gives the upper bound of the range ([2]) . 

Thus, our Universe is created in a hot mixed state, but 
its evolution does not contradict the large-scale structure 
formation. After nucleation from the instanton the Uni- 
verse expands; its radiation dilutes, so that A starts dom- 
inating and generates inflation under an assumption that 
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everywhere above A is a composite field (like an inflaton) 
decaying at the exit by a standard slow-roll scenario. 
The ensemble of universes belongs to a bounded range 
of A = 3H^. Its infrared cutoff is provided by the 
radiation back reaction and survives even in the classical 
limit as _B — > 0. In contrast, the high-energy cutoff at 

A,„ax = 3m|,/2B, ml = 37T/AG, (13) 

is the quantum effect of vacuum energy and the conformal 
anomaly; this generates a new scale in gravity theory. 

We have considered only conformal fields which make 
our model exactly solvable and provide critically im- 
portant positivity of the constant B = 3/3/4, cf. ([7]). 
Moreover, conformal invariance together with the ghost- 
avoidance renormalization renders a particular value of 
the vacuum energy B/2 in ([9|) which yields the upper 
boundary of ^ exactly at the critical point (1/25, B/2) 
of FigO Even if non-conformal fields qualitatively pre- 
serve the whole picture, they are likely to break this re- 
lation. Then if Cq < B/2 all garlands survive, though 
they saturate at A,„ax with a finite temperature. If 
B > Co > B/2, their sequence is truncated at some k. 
Finally, if Co > B the infrared catastrophe occurs again 
- the k = 1 instanton family hits the C-axes at Cq. 

Conformal invariance can be justified as a good ap- 
proximation when conformal particles outnumber non- 
conformal ones. Moreover, their large number N justifies 
a semiclassical expansion by scaling down the range 
Indeed, for a single scalar field the latter is determined 
by Planckian values, Amin ~ 8.99 mp, A^ax ~ 360 mp 
which, however, decrease as 1/A'^ in view of the simple 
scaling C NC, B ^ NB, F{t]o) NF{r]o) and 
ij2 _^ 7j2^jY. Semiclassical expansion can also be justi- 
fied for large i? = 3/3/4 growing with spin, cf.©, because 
the domain ^ with (|13p shrinks to a narrow subplanck- 



ian range when ascending the particle hierarchy. 

Though motivated by the string landscape, all the 
above results hold outside of the string theory context 
and, as a feedback, suggest a long-sought selection mech- 
anism for the plethora of string vacua. Modulo the de- 
tails of a relevant 4D-compactification, this might work 
as follows. For B growing with N and spin, the upper 
scale (jl3l) decreases towards the increasing phenomenol- 
ogy scale, and approaches the latter at the string scale 
TOg where a positive A might be generated by the mech- 
anism of [20|. Our conjecture is that at this scale our 
bootstrap becomes perturbatively consistent, provided 
m?p/B ~ TOg <C TOp, and selects from the string land- 
scape a small subset compatible with observations. 

Our results hold within the Euclidean path integral 
([3]) which automatically excludes Lorentzian config- 
urations possibly existing above the upper boundary 
of pn]) . 4:CH^ > 1. However, one can imagine an 
extended formulation of quantum gravity generalizing 
([3]) to a wider path integration domain. Our conclu- 
sions nevertheless remain true. Indeed the effective 
action scales as Iq ~ —VB, i? ^ 1, and because it 
is negative our landscape at the scale iris is weighted 
by exp{=f^\/B) ~ exp(# mp/ms) S> 1. Therefore it 
strongly dominates over Lorentzian configurations, the 
amplitudes of the latter being 0(1) in view of their pure 
phase nature. Thus, our results look robust against 
possible generalizations of Euclidean quantum gravity. 
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